Nature of well-defined conductance of amine anchored molecular junctions 
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Amine terminated molecules show well behaved conductance in the scanning tunneling micro- 
scope break-junction experimental measurements. We performed density functional theory based 
electron transport calculations to explain the nature of this phenomenon. We find that amines can 
be adsorbed only on apex Au atom, while thiolate group can be attached equally well to under- 
coordinated and clean Au surfaces. Our calculations show that only one adsorption geometry is 
sterically and energetically possible for amine anchored junction whereas three different adsorption 
geometries with very distinct transport properties are almost equally probable for thiolate anchored 
junction. We calculated the conductance as a function of the junction stretching when the molecules 
are pulled by the scanning tunneling microscope tip from the Au electrode. Our calculations show 
that the stretching of the thiolate anchored junction during its formation is accompanied by signifi- 
cant electrode geometry distortion. The amine anchored junctions exhibit very different behavior - 
the electrode remains intact when the scanning tunneling microscope tip stretches the junction. 

PACS numbers; 73.63.-b, 85.65. +h, 72.10.-d 
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I. INTRODUCTION 

Considerable experimental and computational ef- 
forts have been devoted in recent years to molecular 
electronics.^ The main goal in molecular electronics is 
the formation and characterization of nano-electronic de- 
vice in which a molecule plays the role of a conduct- 
ing element. Molecular junctions are formed by plac- 
ing molecules in thin, sub-nanometer gap between metal 
electrodes. It is assumed that under certain conditions 
a single molecule bridges the gap and conducts the elec- 
tricity. Although there has been a significant progress 
towards experimental determination of a single molec- 
ular conductance, there is growing awareness that the 
number of molecules which bridge the gap is not a well 
characterized experimental quantity. Indeed, it is dif- 
ficult to prove experimentally without uncertainty that 
the measured value is the single molecular conductance 
rather than the conductance through molecular array. 

In the experimental technique, which was pioneered 
by Tao at al^t^ and then used by Venkataraman et 
al^^ for amine terminated molecules, the junctions are 
formed by repeatedly crashing an Au scanning tunnel- 
ing microscope (STM) tip into and moving it out of 
Au surface in a molecular solution. The recent exper- 
iment of Venkataraman et al.^ demonstrated that amine 
anchored molecular junctions exhibit well-defined and 
experimentally stable transport properties, namely 1,4- 
benzencdiaminc (BDA) has 2-3 orders smaller variations 
of the conductance between different experimental mea- 
surements than 1,4-benzenedithiolate (BDT)^ This is a 
counterintuitive observation. The thiolate-Au chemical 
bond is considerably stronger than the amine- Au bond. 
Therefore, one may expect that the stronger molecule- 
electrode link would reduce the stochastic switching and 
result into more stable molecular junction. Moreover, 
the weak amine coupling to the electrode implies the 
decreased junction "transparency" for electron current. 
What is the reason for well-defined, stable conductance 



of amine anchored single-molecular junctions? Why do 
amine terminated molecules tend to form single molecu- 
lar junctions? These are the questions which interest us 
here. 

Our paper is organized as follows. In section II, we 
present computational and theoretical methods. In sec- 
tion III, we address three important aspects of the inter- 
face chemistry and their roles in the electron transport: 
(i) possible adsorption sites for amines and thiolates on 
clean and undercoordinated Au(lll) surfaces, (ii) sen- 
sitivity of electron transport to the interface geometry, 
(iii) conductance evolution as the junction is formed and 
stretched by the STM tip. In section III, we also analyze 
and discuss these results and elucidate the nature of the 
well-defined conductance of amine terminated molecular 
junctions. Section IV summarizes the main results of our 
study. 



II. COMPUTATIONAL METHODS 



We used three computer programs in our calculations: 
SIESTA codei to study adsorption of molecules on Au 
surface, ATK package^iS to compute transmission spectra 
for all possible adsorption geometries, and our own code 
for "on-the-fly" conductance calculations to model STM 
break junction experiments. In all our calculations we 
use the same level of electronic structure theory. Namely, 
we apply density functional theory (DFT) with general 
gradient approximation (GGA-PBE) for the exchange- 
correlation potentialiiS, We use numerical single-C with 
polarization (SZP) for Au and the rest of the atoms 
are described by numerical double-^ with polarization 
(DZP) basis sets,! The core electrons are modeled with 
the TrouUier-Martins nonlocal pseudopotentialsf^i 
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A. Conductance calculations for static junctions 

ATK program is based on the combination of the non- 
equilibrium Green's functions (NEGF) and DFTi^i^ We 
describe below the details of the ATK program relevant 
to our applications. The program uses the surface DFT 
calculations to obtain the electrode self-energyi- The ma- 
trix product of the Green's function and the imaginary 
part of the left/right electrode self-energy is used to de- 
fine the spectral density. The non-equilibrium, voltage- 
dependent density matrix is computed from the spectral 
density. Then the density matrix is converted into non- 
equilibrium electron density. The non-equilibrium elec- 
tron density yields the Green's function for the scatter- 
ing region. The Hartree potential is determined by the 
solution of the Poisson equation with voltage-dependent 
boundary conditions. This iterative step [Green's func- 
tion — !■ nonequilibrium electron density — > Green's func- 
tion] is repeated until the self-consistency is achieved. 
Then, the transmission spectrum is calculated by the 
standard equationi^ 



nanowires j 1"'^^ In this model, electrode energy levels are 
broadened to obtain the density of states, which is re- 
lated to the imaginary part of self energy. Our trans- 
port implementation avoids the computationally expen- 
sive Cauchy principal value numerical integration and 
uses the Lorentz broadening instead, as was suggested 
in our non-orthogonal Wannier-type atomic orbital elec- 
tron transport calculations The Lorentzian broadening 
cr in the electrode density of state is equivalent to positive 
infinitesimal a in the electrode Green's function>i^ 

[(ii; + *^)s-H]-i , (1) 

where S is the overlap matrix and H is the electrode 
Hamiltonian. The conductance of monoatomic gold 
chain attached to two gold clusters is Gq — 2e^/h^ 
and the broadening parameter a is fitted to reproduce 
this value. This fitting leads to a—O.b eV and we use 
this value in all our calculations. This computational 
method enables us to perform conductance calculations 
for hundreds of sample geometries from a junction evo- 
lution trajectory. 



B. Junction formation and stretching: on-the-fly 
conductance calculations 

We would like to model the experiments where the 
molecular junctions are formed by repeatedly crashing 
an Au STM tip into and pulling it out of Au sur- 
face in a molecular solutionj^i^i^i^i^ To model these ex- 
perimental conditions we developed the computational 
method which calculates molecular conductance " on-the- 
fly" . Namely, we performed conductance calculations 
along the trajectory, which is computationally obtained 
by pulling the molecule from the Au electrode. 

The geometry evolutions of the molecular junctions are 
obtained by constrained geometry relaxation methodJ^^ 
The molecular junction is aligned along the z-axis. The 
z coordinate of the top C atom is incremented by 0.05 
A on every step and this stretch of the junction is fol- 
lowed by the constrained geometry relaxation. In the 
constrained geometry relaxation, the z coordinate of the 
central bottom Au atom of the electrode cluster and the 
top C atom of the phenylaime or phenylthiolate molecule 
are fixed and the rest of the system is fully optimized. 
The upper electrode is the mirror image of the bottom 
one. The junction evolution is represented by the trajec- 
tory composed of the optimized geometries. It becomes 
the true molecular dynamics trajectory in the limit of ex- 
tremely slow pulling speed. Previous studies found that 
the constrained geometry relaxation method gives the re- 
sults consistent with Car-Parrinello molecular dynamics 
simulation of pulling molecules from surface. ^^^^ 

For every step along the trajectory the conductance 
is computed by the method, which we implemented 
in the SIESTA program.^ We compute the self en- 
ergy by the level broadening model, which has been 
successfully applied to molecular junctions and metal 



III. RESULTS AND DISCUSSION 

A. Phenylamine and phenylthiolate adsorption on 
clean and undercoordinated gold surfaces 

First principles studies of adsorption of different amine 
or thiolate terminated molecules on Au(lll) surfaces 
have been reported. "'^^i-'^^'^^i^^ It was predicted that am- 
monia is adsorbed on the atop siteJ^ The cluster model 
calculations show that the adsorption energy depends 
significantly on the employed theoretical method^ 
Phenylthiolate is attached to the bridge site and the 
phenyl ring is strongly tilted with respect to the Au 
surface/^ Although previous studies are very informa- 
tive, a comparative study of adsorption of amines and 
thiolates on Au surfaces in the context of electron trans- 
port calculations is necessary. 

We represent Au(lll) surface by unit cell, which is 
periodically repeated in two directions. The unit cell 
consists of five (4 x 4) layers of Au atoms. The under- 
coordinated gold surface is modeled by adding an extra 
apex atom above the Au(lll) hollow site. 

We perform geometry optimizations for phenylamine 
and phenylthiolate molecules adsorbed on clean and 
undercoordinated Au(lll) surfaces. A final optimized 
structure is sensitive to initial guess for geometry op- 
timization. We choose various starting geometries to 
explore all possible adsorption sites. FiglT] shows atop, 
bridge and hollow sites on Au(lll) surface. There are 
two different hollow sites on Au(lll) surface: fcc-hoUow 
and hcp-hollow sites. We initially align the phenyl ring 
of the phenylamine and phenylthiolate along the Au-Au 
or the Au-hoUow directions. The phenyl ring of pheny- 
lamine and phenylthiolate can be perpendicular or tilted 
with respect to the Au surface. For tilted phenylamine. 



FIG. 1; Diagram which shows possible adsorption sites on 
Au(lll) surface. The Au atoms are represented by circles of 
different sizes. The smaller the spheres, the farther from the 
layer the Au atoms are. 



TABLE I: Adsorption geometries of phenylamine and 
phenylthiolate. A'au is the number of the nearest neighbor 
Au atoms for the N atom or the S atom. cJau-n and dAu-s 
are the mean distances for the A^au nearest neighbors of the 
N atom and the S atom respectively. 6 is the angle between 
the phenyl ring and the surface (Fig[T]a). 

phenylamine 

adsorption geometry A'au dAu-N(A) S(degree) AE (eV) 



FigEl(b) 
FigEKc) 
Figll(d) 
FigEKe) 



1 


2.90 


0.83 


1.92 


2 


3.05 


87.5 


0.68 


1 


2.34 


87.9 


1.18 


1 


2.30 


12.0 


0.92 



phenylthiolate 
adsorption geometry Aau dAu-s(A) 6'(degree) AE (eV) 

Figl2](f) 
Figl2](g) 
FigEl(h) 
Figl2](i) 
Fig|2](j) 



2 


2.56 


8.12 


3.64 


3 


2.44 


89.2 


2.57 


3 


2.51 


63.7 


2.72 


1 


2.31 


31.5 


2.72 


1 


2.31 


51.3 


2.75 





FIG. 2: adsorption geometries for phenylamine (b-e) and 
phenylthiolate (f-j) on the clean (b,c,f,g,h) and undercoordi- 
nated (d,e,i,j) Au(lll) surfaces, (a) defines the angle 6, which 
is used in Table I. 



the initial orientation is chosen in such a way that the 
N lone pair points perpendicular to the surface. We 
choose angle 30° for the initial phenylthiolate tilted ori- 
entation. The coordinates of the Au atoms in the two 
bottom layers are constrained to theoretical bulk geom- 
etry with lattice vector 4.132 A . This lattice vector is 
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obtained by performing a series of bulk calculations with 
different lattice parameters and fitting the energies to 
the third-order Birch-Murnaghan equation of state.— We 
performed two subsequent geometry optimizations to get 
equilibrium structures for the adsorbed molecules. In the 
first geometry optimization, the N atom or the S atom 
(and the Au apex atom when the surface is undercoor- 
dinated) are fixed on a specific adsorption site and can 
only move in the z direction, which is perpendicular to 
the surface. In the second one, all atoms except the two 
bottom Au layers are fully relaxed. 

Having obtained the optimized geometries for the 
molecules adsorbed on the surface, we compute adsorp- 
tion energy (AE) by subtracting the energy of the total 
system from the sum of energies of the surface and the 
molecule.— Figl2] and Table I summarize the results of 
the geometry optimizations. Due to the failure of DFT 
GGA exchange-correlation functionals to account for van 
der Waals interactions, the actual AEs may be somewhat 
higher than the calculated values. AE is mainly deter- 
mined by the orientation of the phenyl ring with respect 
to the surface. When the phenylamine is adsorbed on 
the clean Au (111) surface its phenyl ring either lies al- 
most flat on the surface (geometry (b) on Fig. [2|) or 
stays vertical to the surface (geometry (c) on Fig. 
AE of geometry (b) is about 1.2 eV larger than that of 
geometry (c). This can be expected since there is strong 
attractive interaction between the molecular 7r-electrons 
and the metal surface image charges, which favors par- 
allel to the surface phenyl ring orientation. The phenyl 
ring lies above the three- Au triangle on the Au(lll) sur- 
face. The differences in the adsorption positions of the 
N atom cause only small AE variations of about 0.1 eV, 
if the phenylamine is attached vertically to the surface. 
When the phenylamine is adsorbed on undercoordinated 
Au(lll) surface, the phenylamine adopts vertical (geom- 
etry (d) on Fig. m or tilted (geometry (e) on Fig. ^ 
configurations. The N lone-pair points to the apex Au 
atom in both cases. AE of the vertical phenylamine (1.18 
eV) is 0.26 eV larger than that of the tilted one, and it 
is 0.50 eV larger than AE of the most stable non-lying 
adsorption geometry on clean Au(lll) surface. 

The geometry optimizations of phenylthiolate on clean 
Au(lll) surface also yield adsorption structures where 
the molecule either lies almost flat on the surface (f), or 
it is vertical (g), or it is tilted (h) with respect to the 
surface. The lying phenylthiolate is not perfectly paral- 
lel to the surface, because the attraction between the S 
atom and the surface pulls the atom closer to the sur- 
face. In agreement with previous calculations,^" the S 
atom appears near the bridge site in the lying geometry. 
If phenylthiolate is adsorbed on the apex Au atom, the 
molecule is always tilted with respect to the surface. For 
phenylthiolate, there are two possible tilted geometries. 
The phenyl ring plane is tilted in the first geometry ((i) 
on Fig. [21), and the second geometry has the vertical 
phenyl ring but the tilted S-C-C axis ((j) on Fig. 
The AE difference of these two possible tilted geometries 
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FIG. 3: Transmission spectra of BDA (upper panel) and BDT 
(lower panel) computed for different adsorption geometries. 
Fermi energy of the electrodes is shifted to zero. The junction 
geometries (d), (h),(i) and (j) refer to Table I and Fig[2] 

is within 0.1 eV. The lying adsorbed molecule also gives 
the largest AE for phenylthiolate. But unlike pheny- 
lamine, the tilted adsorption geometries of phenylthiolate 
on clean and undercoordinated Au(lll) surface have very 
similar AEs, which are about 0.9 eV smaller than that 
for the lying geometry and 0.15 eV larger than AE for 
the vertical geometry. 



B. Electron transmission dependence on 
adsorption site geometry 

Figl21 and Table I list the most probable geometries 
for the adsorption of phenylamine and phenylthiolate on 
Au(lll) surface. These configurations are obtained from 
the geometry optimizations described in the previous sec- 
tion. The lying on the surface molecular junctions (Fig|5] 
(b), (e), (f)) are excluded from the transport calculations, 
since it is difficult to form them experimentally. Our cal- 
culations show that their conductance is close to or even 
larger than the conductance of a gold monoatomic wire. 
Go = 2e^//i. It means that the lying on the surface junc- 
tions are also undetectable in the STM break-junction ex- 
periments. Three different geometrical arrangements (h) , 
(i), (j) are energetically and sterically possible for BDT 
junction to the Au surface, whereas only one geometry 
(d) is available for the BDA junction. Given these opti- 
mized adsorption geometries, the molecular junctions are 
constructed by mirror addition of the other electrode and 
the corresponding anchoring group. The conductances of 
these four junctions are calculated by NEGF/DFT as is 
implemented in the ATK package. ^'^ Transmission spec- 
tra of the BDA and the BDT junctions are plotted in Fig. 
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Stretching (A) 



FIG. 4: a) Snapshots of the junction geometries, b) Conductance traces of the BDA (left curve) and the BDT (right curve) 
junctions. The starting point of conductance traces for the BDT junctions is shifted along the x axis to 7.5 A . Inset; the 
starting junction geometries. 



[3l The BDA transmission is featureless in the vicinity of 
the Fermi energy. Fig[3] shows very distinct transmission 
spectra for three geometrically different BDT jimctions 
which have similar AEs. The BDT junction with geom- 
etry (h) does not have any transmission resonances near 
the electrode Fermi energy. The junction with geometry 
(i) shows two resonances: one broad resonance is below 
and one resonance is above the Fermi energy. The most 
"transparent" for electrons junction is the junction with 
geometry (j). It has one peak located at the Fermi en- 
ergy. The transmission spectra for the BDT junctions 
are in quantitative agreement with the results of density 
functional theory based electron transport calculations 
reported by other authors i^^i^^'^^ 
The Landauer formula3^ 

G = G^T{Ef) (2) 
relates the conductance G with the transmission func- 



tion T{Ef) computed at the Fermi energy Ep of the 
electrodes. This linear relation is valid when the applied 
voltage bias is small.— It is sufficient for our applications 
since very low voltage (10-25 mV) is used in STM break 
junction experimentsi^i^. Therefore the conductance for 
BDA and BDT junctions can be extracted from the trans- 
mission spectra plotted on Fig. [S] The conductance of 
the BDA junction is 0.038 Gq. The experimentally mea- 
sured conductance for BDA junction is 0.0064 Gq (re- 
ported in Refi^) that is 6 times smaller than the com- 
puted value. The reason for this discrepancy is not clear. 
The BDT junction shows a wide range of conductance 
values: 0.1 Go for geometry (h), 0.2 Go for geometry (i), 
and 0.7 Go for geometry (j). Such variations are qualita- 
tively consistent with the experimental observation that 
the mean conductance variation for amine anchored junc- 
tion is much smaller than that for thiolate linked junction 
in different experimental conductance tracesi^ 
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Our study demonstrates that BDA adsorbs on apex 
Au atom, while BDT can be adsorbed on both clean Au 
surface and apex Au atom. It means that there are more 
geometrical constraints for BDA than for BDT in form- 
ing junctions. In addition, the N lone-pair leads to the 
stronger orientation preference for the amine-Au bond 
than that for the thiolate-Au linkage. Therefore, it is 
sterically more difficult to get multiple BDA molecules 
in a junction. Based on this picture, we suggest that 
the recently studied dithiocarbamate and dithiocarboxy- 
late anchoring groups which are even more con- 
strained orientationally, would provide stable and well- 
defined junctions for experimental measurements. 



C. Conductance Trace Calculations 

We model STM break-junction experiments by our 
program for "on-the-fly" conductance calculations. We 
use the most stable non-lying adsorption geometries of 
phenylamine and phenylthiolate ((d) and (j) respectively 
on Figl2) to prepare the cluster models for the conduc- 
tance trace calculations. The atomic coordinates of the 
10- atom gold cluster and the molecule are extracted from 
the optimized adsorption geometry. Then we perform 
150 constrained geometry relaxations for each cluster- 
molecule system as is described in section II B. This se- 
ries of the constrained geometry relaxations yields the 
trajectories which describe the evolution of the junction 
when the molecule is pulled from the Au surface by the 
STM tip. Each trajectory uniformly spans the junction 
evolution from the equilibrium stretch A to 150x0.05 
A= 7.5 A. FigHfa) shows representative snapshots of the 
geometry evolutions for the BDA and BDT junctions. 

The junctions undergo completely different geometry 
evolutions as they are stretched. The pulling of BDT 
causes the significant distortion of the Au cluster, which 
represents the electrode. At the same time, the geome- 
try of the Au cluster is not significantly affected by the 
stretching of the BDA junction. This can be understood 
by comparing amine-Au and thiolate-Au bond strengths. 
The amine-Au is a weak bond. When BDA is pulled from 
the Au electrode, the bond between the amine and the 
apex Au atom ruptures before the molecule is able to 
involve the other Au atoms from the cluster into its mo- 
tion. The interaction between thiolate and Au is consid- 
erably stronger than Au-Au interaction, and it involves 
all Au atoms into the motion as the BDT junction is 
stretched. Therefore, the electrode geometry evolution 
for the BDA junction is characterized by two states (be- 
fore and after the amine-Au bond rupture), whereas the 
BDT junction has different electrode geometries for every 
stretching step. 

Junction conductance traces and the geometry evolu- 
tion trajectories are calculated simultaneously (see sec- 
tion II B for detail). Fig. [5] shows the conductance 
traces for the both jimctions. The conductance of BDA 
only slightly fluctuates around several hundredth of Gq 



at early stage of the pulling, then it drops abruptly when 
the amine-Au linkage ruptures. Oppositely, the BDT 
junction undergoes almost two orders of magnitude vari- 
ation in the conductance before the molecule is detached 
from the electrode. This difference of the behavior of 
the conductance variations is consistent with the exper- 
imental observations.^ We also noticed that the conduc- 
tance of the BDT junction is increased as the junction 
is stretched by the STM tip. This behavior of the the 
conductance trace is in agreement with Xue and Ratner 
observation that the the conductance of BDT may grow 
with the increase in the molecule-electrode separation.— 
It is also consistent with Ke, Baranger and Yang calcula- 
tions which show that the presence of an additional Au 
atom at each of the two contacts will increase the BDT 
conduct ancei^ 



IV. CONCLUSION 

We performed a computational study to elucidate the 
nature of the well-deflned conductance of amine termi- 
nated molecular junctions. We analyzed possible ad- 
sorption geometries of amines and thiolates on clean and 
undercoordinated Au(lll) surfaces. Electron transport 
properties of BDT and BDA molecules are computed for 
static junctions as well as "on-the-fly" to model STM 
break junction experiments. The main results of our 
study can be summarized as follows: 

1. Amines adsorb only on undercoordinated Au sur- 
face whereas thiolates adsorb equally well on clean 
and undercoordinated Au surfaces. AE of the verti- 
cal phenylamine on undercoordinated Au(lll) sur- 
face is 0.50 eV larger than AE of the most stable 
non-lying adsorption geometry on clean Au(lll) 
surface. AEs of phenylthiolate on on clean and un- 
dercoordinated Au(lll) surface are very similar (~ 
2.7 eV). 

2. The number of accessible geometries for BDA junc- 
tion on Au surface is much smaller than that for the 
BDT. Only one geometry is sterically /energetically 
accessible for BDA junctions whereas three differ- 
ent geometries are possible for BDT. It suggests it 
is more likely for the amine terminated molecules 
to form single molecule junctions than for thiolate 
terminated molecules. 

3. Conductance of static BDA junction is 0.0064 Gq 
whereas conductance of BDT varies from 0.1 to 0.7 
Go depending on the adsorption geometry. 

4. The conductance of BDA only slightly fluctuates 
at the early stage of the junction formation, then 
it drops abruptly as the amine-Au linkage rup- 
tures. Oppositely, the BDT junction undergoes al- 
most two order of magnitude variation in the con- 
ductance before the molecule is detached from the 
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electrode. This is because the electrode geometry 
is significantly distorted by the thiolate anchored 
junction stretch and is not affected by the stretch 
of the amine anchored junction. 
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